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Red variables in the OGLE-II data base - III. Constraints on the 
three-dimensional structures of the LMC and SMC 
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ABSTRACT 

We present an analysis of the 3-D structure of the Magellanic Clouds, using period- 
luminosity (P-L) relations of pulsating red giants in the OGLE-II sample. By interpreting 
deviations from the mean P-L relations as distance modulus variations, we examine the three- 
dimensional distributions of the sample. The results for the Large Magellanic Cloud, based 
solely on stars below the tip of the Red Giant Branch, confirm previous results on the inclined 
and possibly warped bar of the LMC. The depth variation across the OGLE-II field is about 2.4 
kpc, interpreted as the distance range of a thin but inclined structure. The inclination angle is 
about 29°. A comparison with OGLE-II red clump distances revealed intriguing differences 
that seem to be connected to the red clump reddening correction. A spatially variable red 
clump population in the LMC can explain the deviations, which may have a broader impact 
on our understanding of the LMC formation history. For the Small Magellanic Cloud, we 
find a complex structure showing patchy distribution scattered within 3.2 kpc of the mean. 
However, the larger range of the overall depth on every line-of-sight is likely to smooth out 
significantly the real variations. 
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1 INTRODUCTION 

The close proximity of the Magellanic Clouds offers an impor- 
tant opportunity to study in detail the dynamics and composition 
of these galaxies. One of the important issues that affects our in- 
terpretation of the observable parameters is the three-dimensional 
structure of the Clouds. The Large Magellanic Cloud (LMC) has 
long been regarded as a thin flat disk seen nearly face-on, but with 
the east side being closer than the west (Westerlund 1997; see also 
van der Marel 2004 for an updated review). Recent determinations 
of the tilt angle place it in the range 30-35° (van der Marel & Cioni 
2001, Olsen & Salyk 2002, Nikolaev et al. 2004). However, one can 
find a significantly broader range in the literature (25^15°, Wester- 
lund 1997), which can at least partly be explained by nonplanar 
geometry of the inner LMC (Nikolaev et al. 2004). In contrast, the 
structure of the Small Magellanic Cloud (SMC) has proved difficult 
to define, mostly because of its large inclination angle. Different 
distance indicators implied a line-of-sight depth from 7 kpc (O, B 
and A type stars; Azzopardi 1982) to 6-12 kpc (star clusters; Crowl 
et al. 2001) and 20 kpc (Cepheids; Mathewson et al. 1986, 1988). 



* E-mail: l.kiss@physics.usyd.edu.au 

f On leave from University of Szeged, Hungary 



It is well established that the asymmetric appearence of the SMC 
is due almost exclusively to the young population of stars, since 
the older stellar population shows a very regular distribution (e.g. 
Zaritsky et al. 2000, Cioni et al. 2000, Maragoudaki et al. 2001). 

Time-series observations from microlensing projects such as 
MACHO (Alcock et al. 2000) and OGLE (Udalski et al. 1992), 
combined with data from two recent near-infrared surveys, 2MASS 
(Skrutskie 1998) and DENIS (Epchtein 1998), have provided an 
unprecedented new wave of high-quality data and a dramatic im- 
provement in statistical analyses. The order-of-magnitude increase 
in the numbers of well-measured stellar distance indicators has 
allowed detailed searches for small distance modulus variations 
across the face of the Clouds. Most recent examples for the LMC 
include the Cepheid period-luminosity (P-L) relations (Nikolaev 
et al. 2004), the Asymptotic Giant Branch (AGB) luminosity dis- 
tribution (van der Marel & Cioni 2001) and the luminosity of core 
helium-burning red clump stars (Subramaniam 2003, 2004). The 
basic idea was the same in all these studies: the observed ap- 
parent magnitude deviations were translated to distance modulus 
variations. The analysed fields ranged from the central 4.5 square 
degrees of the LMC bar (Subramaniam 2003, 2004) to over 300 
square degrees, covering the outest parts of the disk (van der Marel 
& Cioni 2001). 
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Table 1. The inverse regression coefficients of the fitted P-L relations. They 
have the form log P = a 1 X Kg + b' . 'Number' refers to the number 
of stars used in the line fit for that P-L relation. 



P-L relation" 


a' 


b< 




Number 




LMC 








Rs (A") 


-0.237 ±0.002 


4.340 ± 0. 


.029 


2642 


R 2 (B-) 


-0.269 ± 0.004 


4.931 ±0. 


.045 


1634 




SMC 








Rs(A-) 


-0.206 ± 0.006 


4.030 ± 0. 


.120 


229 


R 2 (B") 


-0.214 ±0.006 


4.330 ± 0. 


.130 


117 


30 (A+) 


-0.206 ± 0.006 


4.038 ± 0. 


.070 


133 


20 (B+) 


-0.214 ±0.006 


4.349 ± 0. 


.069 


218 


lO(C') 


-0.238 ±0.006 


4.821 ±0. 


.069 


260 


F(C) 


-0.222 ±0.007 


4.932 ± 0. 


.078 


405 


L 2 (D) b 


-0.170 ±0.005 


4.799 ± 0. 


.062 


534 


L 2 (D) c 


-0.170 ±0.005 


4.810 ±0. 


.160 


105 



a - abbreviations by Ita et al. (2004a) in parentheses 
b - above the tip of the Red Giant Branch 
c - below the tip of the Red Giant Branch 



In this Letter we present the first attempt to use red giant 
P-L relations to constrain the three-dimensional structures of the 
Clouds. Since the discovery of multiple P-L relations of long- 
period variables in the LMC (Wood et al. 1999, Wood 2000), there 
has been strong interest in pulsating red giants and many indepen- 
dent analyses have been published (Cioni et al. 2001, 2003; Noda 
et al. 2002, 2004; Lebzelter, Schultheis & Melchior 2002; Kiss & 
Bedding 2003, 2004 - Paper I and Paper II; Ita et al. 2004ab; Groe- 
newegen 2004; Schultheis, Glass & Cioni 2004; Soszyhski et al. 
2004). Interestingly, none of these studies considered the geometric 
effects on the apparent magnitudes, as has been done, for instance, 
by Wray, Eyer & Paczyhski (2004), who analysed short-period red 
giants in the galactic Bulge. The main aim of the present paper is 
to address this issue and draw conclusions on the line-of-sight dis- 
tance variations of the Clouds, as traced by pulsating red giants. 



2 DATA ANALYSIS 

In this work we used OGLE-II periods and 2MASS K s magnitudes 
of red giant variables in the LMC and SMC, as discussed in Papers 
I and II, where all relevant details on data reduction can be found. 
Our key assumption here is that the vertical scatter about any par- 
ticular P-L relation contains information on the relative distances 
to individual stars. This allows us to examine the magnitude dif- 
ferences from mean P-L relations as a function of position on the 
sky. Although this sounds simple, we have to ask to what extent 
this assumption is valid. There are three aspects to the question: 
( i) star-to-star extinction variations due to interstellar dust; ( ii) the 
contribution of stellar variability to the vertical scatter in the P-L 
relations; (Hi) the intrinsic width of the P-L relations. 

(i) The influence of dust absorption in the near-infrared was 
discussed by van der Marel & Cioni (2001) for the LMC, who con- 
cluded that it was only a few hundredths of a magnitude in Ks 
and could be neglected. This has later been criticised by Nikolaev 
et al. (2004), who preferred the multi-wavelength approach that in- 
volved both distance modulus and extinction variations. Their fig. 
7 shows that most of the individual Cepheid reddenings scattered 
within 8E(B — V) ~ ±0^1, which can be translated to a random 



error 5K S ~ ±0™03, using the extinction law in Schlegel et al. 
(1998). However, there are also two facts which have to be con- 
sidered: 1. the Nikolaev et al. sample covered a much larger field 
of view, where larger extinction variations can occur; 2. Cepheids, 
as young supergiant stars, are often found close to star forming re- 
gions, where the local dust amount might be much higher than the 
average. Individual reddenings have also been determined by Sub- 
ramaniam (2003), who analysed exactly the same OGLE-II field of 
view as us (Udalski, Kubiak & Szymanski 1997), allowing a direct 
comparison. She found A ma xE(V — I) « 0? 1 029 as the maximum 
reddening dispersion along the bar of the LMC, which corresponds 
to 8K S ~ ±0^01 extinction dispersion in Ks- From these num- 
bers, we conclude that for each star, the uncertainty in Ks due to 
extinction is likely to be less than 0^03, possibly around 0^01, 
which is indeed negligible if we can take averages over hundreds 
of stars. 

(ii) Since 2MASS magnitudes are single-epoch measure- 
ments, intrinsic variability can introduce some uncertainty. For that 
reason, we decided to use only small-amplitude first ascent Red 
Giant Branch (RGB) stars in the LMC, of which several thousands 
were identified in Paper I. Their amplitudes in / range from 0^005 
to 0^02, so that 2MASS Ks measurements should be within 
0? 1 01-0? 1 02 of the mean values. The situation is less favourable in 
the SMC, for which the smaller number of stars (about 3,200 were 
analysed in Paper II) forced us to use most P-L relations, including 
both RGB and AGB variables. The latter may have amplitudes up 
to several tenths of a mag in the infrared (Whitelock, Marang & 
Feast 2000). 

( Hi) The horizontal width of the P-L relations is affected by er- 
rors arising from period determination and instrinsic scatter caused 
by astrophysics (e.g., mixing different populations of stars or ran- 
dom and time-dependent excitation of pulsations). Considering P- 
L ridges R 2 and R3, we compared the relations presented in Paper I 
with those determined by Wood (2000) and Soszynsky et al. (2004), 
both based on about 8 years of observations (twice as long as 
OGLE-II). If the widths were dominated by errors in period deter- 
mination caused by the short time-span of data, the longer datasets 
should have led to significantly tighter P-L relations. Somewhat 
surprisingly, this is not the case: we found the same ~0. 1 dex width 
in log P in all three studies, which implies that the natural width of 
the relations is probably resolved by OGLE-II alone. We note that 
errors in period determination can also arise from the fact that the 
majority of stars show multiply periodic behaviour, so that period 
uncertainty does not necessarily scale inverse-proportionally with 
the full time-span. In this work we assumed that the scatter in the 
P-L relations is partly due to a random error in period, which will 
be averaged out in a large sample, and that there is a vertical scatter 
due to a statistical error in the distance moduli estimates. 

To calculate distance modulus variations, each star in our sam- 
ple was classified as lying on one of the P-L relations. To minimize 
the effects of spurious periods, we only used the dominant period 
for each star. For the classification, we drew straight line bound- 
aries between the different relations, similarly to Ita et al. (2004b). 
The results do not depend heavily on the exact definition of these 
lines because the majority of the stars lie away from the boundaries. 
To distinguish between AGB and RGB stars, we adopted TRGB 
magnitudes from Paper II. We then made linear least-squares fits 
to each of the P-L relations. Because of the non-zero width and 
the rhomboid-shaped P-L ridges, these fits were performed as in- 
verse regressions in the form of log P — a' x Ks + b' (Table 
0, which were then converted back to usual P-L relations in form 
of Ks — ax log P + b. For a star at a given period, the verti- 
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Figure 1. The binned K$ magnitude differences for various Right Ascen- 
sions, compared with red clump results by Subramaniam (2003). Larger 
values correspond to those parts of the LMC that are further from us. 

cal difference between its observed Ks magnitude and the linear 
fit was taken as the distance modulus of that star relative to the av- 
erage distance of the host galaxy. Because of the relatively large 
scatter in these distance moduli, caused by the issues discussed 
above, it was helpful to bin the distance data to see the underlying 
trends. When combining the results from different P-L relations, 
each measured distance modulus was weighted by the overall scat- 
ter (combined with the number of points) in its P-L relation. This 
gave more weight to the distances based on the P-L relations that 
had the tightest overall correlations. 



3 RESULTS 

3.1 Structure of the LMC 

The results presented here for the LMC used only the R3 and 
R2 period-luminosity relations from Paper I, i.e. those below the 
TRGB (sequences A~ and B~ in Ita et al. 2004a). They had by 
far the tightest relations and gave the best distance measurements; 
4,276 stars were used for the LMC. 

The OGLE-II data are limited to the region around the bar 
of the LMC, which lies approximately in the east-west direction. 
Therefore, the obvious way to bin the data was in Right Ascen- 
sion. This is the same approach as used by Subramaniam (2003, 
2004), who analysed dereddened OGLE-II red-clump magnitude 
variations in terms of relative distance variations. She has kindly 
provided her data and we compare these with our results in Fig. 
Q To overlay the Subramaniam Io distance moduli, we subtracted 
their weighted mean value (((Io) = 18™ 16). Larger values in Fig. 
Q refer to those parts of the LMC that are further from us. The 
uncertainty in each RA bin was calculated from the scatter of the 
distance moduli that were combined to make that point. 

It is clear from Fig. Q that red giant P-L relations alone re- 
veal that the bar is inclined away from us, running east-west. This 
is in perfect agreement with the general view of the LMC (West- 




Figure 2. A 3-D representation of the LMC. The lighter regions are closer 
to us and the darker regions are further away. 

erlund 1997). The magnitude variations can be translated to dis- 
tance variations by assuming a distance modulus to the LMC of 
(m - M) = 18™5 (Alves 2004a), corresponding to 50.1 kpc 
mean distance. The maximum distance modulus range we observe 
in the LMC is Of l±0f03 (calculated from a linear fit between 
RA=77° and 87°), which can be translated to 2.4±0.7 kpc distance 
variation. We interpret this as the distance range of a thin but in- 
clined structure, where the thickness of the LMC bar is assumed 
to vary less than the mean distance along any given line-of-sight. 
Using Fig.Qand the assumed LMC distance, an inclination angle 
of 29 degrees (with a formal error of several degrees) can be de- 
termined. However, note that OGLE-II data are not well suited for 
inclination angle determination because the field of view is highly 
linear so the position angle of the node, along which the inclination 
angle has to be measured, is weakly constrained. In this sense, the 
inclination angle we determined is a lower limit to the real value. 
Nevertheless, our 29° is in good agreement with recent inclination 
angle determinations (e.g., van der Marel & Cioni 2001, Nikolaev 
et al. 2004). 

There is also evidence for substructures within the bar as de- 
viations from a straight line in Fig. Q Similar substructures were 
recently found by Subramaniam (2003, 2004), who interpreted red 
clump magnitude variations as evidence of a misaligned secondary 
bar within the primary one. There appears to be good consistency 
between the two data sets, with the distance values mostly agree- 
ing to within la. However, there are also some discrepancies, most 
prominently at RA~ 79° and RA~ 83 - 85°. At these position 
there are two ~0^05 dips in the Subramaniam data but not in our 
data (although we have fewer data points, so that the 2-3 a differ- 
ences may arise purely from statistical fluctuation). There are two 
reasons to believe that our data in these regions are to be preferred. 
Firstly, optical and narrow-band Ha images show prominent HII 
regions in the quoted positions (e.g. giant shells No. 51, 54 and 60 
at RA~ 79° and No. 77 at RA~ 84° in Kim et al. 1999), which 
suggests that the difference may indicate undercorrected extinction 
in the J-band data used by Subramaniam (2003). Secondly, the two 
RAs where we found the largest deviations coincide exactly with 
the positions where the average E(V — /) reddening values show 
sudden jumps in fig. 2 of Subramaniam (2003). Both facts imply 
that the red clump method (developed by Olsen & Salyk 2002 and 
used by Subramaniam 2003) may be affected by a systematic error 
that is related to the red clump reddening determination. 
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Figure 3. The binned K$ magnitude differences for various Right Ascen- 
sions. The larger values refer to those parts of the SMC that are further from 
us. 



The red clump method assumes the red clump has constant 
(V — /)o colour everywhere in the LMC. Since we used reddening- 
insensitive if-band data for red giants, the discrepancy with clump 
results implies the intrinsic (V — I)o colour of the clump must be 
redder than the assumed constant value ((V — J)o=0 1 "92, Olsen 
& Samolyk 2002), and this probably shows a change in the mean 
red clump population. Alves (2004a) showed that population effect 
on model red-clump colour-magnitude diagrams can significantly 
affect distances derived from V and / band red clump data. Our 
findings can be explained by a spatially varying clump population 
in the LMC, which may have broad implications for understanding 
the LMC bulge/halo/disk structure and formation history, as dis- 
cussed recently by Alves (2004b) and Zaritsky (2004). 

To search for structure in the Declination direction we cre- 
ated a 3-D representation of the LMC in which the relative distance 
modulus is shown in grey scale (Fig.|2j. This figure was made us- 
ing a two-dimensional averaging with a Gaussian weight-function 
(FWHM — 8.'6). This representation makes the depression at 
RA~ 77° quite prominent. The overall appearance shows the tilt 
of the bar, with no measurable systematic trend along Declination 
within the OGLE-II field of view. 



Figure 4. A 3-D representation of the SMC. The lighter regions are closer 
to us and the darker regions are further away. 



larger depth of the SMC (cf. Sect. 1). The latter effect is particularly 
prominent in fig. 9 of Ita et al. (2004a): all ridges, without excep- 
tion, are significantly thicker, including those of the fundamental 
and first overtone Cepheids. For that reason it would be misleading 
to conclude from Fig. [3] that no relative distance modulus varia- 
tions occur in the SMC. There are hints of substructures but their 
interpretation is difficult because the mean distance modulus is a 
complicated function of stellar distribution along the line of sight. 
Formally, the given distance modulus variations correspond to a 
distance range of 3.2±1.6 kpc, adopting (m - M)o = 18™94 (Pa- 
per II), which, contrary to the LMC, is likely to be more affected 
by the overall depth on every line-of-sight. A better depth measure- 
ment would require a study of the vertical scatter along individual 
P-L ridges, but for that purpose, these red giant P-L relations in 
the SMC are less useful because they are too close to each other in 
the P-L plane, which makes it more difficult to distinguish between 
members of different ridges. 

Finally, the 3-D representation of the SMC (Fig.[4j shows that 
both the north-eastern and the south-western corners of the OGLE- 
II field are slightly closer to us, while there is a system of depres- 
sions along the northern boundary. Interestingly, the two deepest 
"holes" (at RA~14° and 11°) coincide exactly with the two ma- 
jor concentrations of red giant (RGB and AGB) stars in figs. 6-7 in 
Cioni et al. (2000), thus they represent significant extensions of the 
main body of the SMC, inclined further away from us. 



3.2 Structure of the SMC 

In the SMC the smaller number of stars led us to use more P-L 
relations to reduce the errors. We included P-L ridges R3, R2, 30, 
20, 10, F and L 2 (A - , B" , A+ , B+ , C', C and D in Ita et al. 2004a); 
these were the only relations for which reasonable boundaries could 
be defined, due to the larger scatter in the SMC. As for the LMC, 
we weighted the results by the corresponding scatter in each P-L 
relation, so that better-defined relations received higher weight. 

The resulting distance modulus variations are plotted in Fig.|5| 
The uncertainties are considerably larger than for the LMC, which 
is partly due to having fewer stars and partly due to the known 



4 CONCLUSIONS 

We have demonstrated the usefulness of pulsating red giants as dis- 
tance indicators in external galaxies by measuring the 3-D struc- 
tures of the Magellanic Clouds. In particular, the variables below 
the tip of the Red Giant Branch are numerous and their narrow P- 
L relations are real competitors to those of the Cepheid variables. 
About a year of continuous monitoring gives very well-defined P- 
L relations for those stars, which have periods between 15 and 50 
days. Because of the much lower amplitudes, intrinsic variability 
introduces negligible vertical scatter in the relations. This is, of 
course, a major observational obstacle: one has to reach the several 
millimagnitude precision over the whole photometric monitoring 
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period, which should be as continuous as possible. For that rea- 
son, red giant P-L relations in general can be considered as giving 
supplementary information to other distance indicators, although 
high-amplitude red giants have already been used for measuring a 
stand-alone extragalactic distance to NGC 5128=Cen A (Rejkuba 
2004). 

From an analysis of almost 4,300 RGB variables in the LMC 
we determined spatial depth of the LMC bar as a function of ce- 
lestial position. The results are in good agreement with those based 
on other distance indicators (Cepheids, red clump stars, AGB lu- 
minosity distributions), but are perhaps less affected by dust. We 
also found possible evidence for spatial variations in the red clump 
population in the LMC, which suggests that the assumption of con- 
stant (V — I)o red-clump colour over the whole LMC may not be 
valid. For the Small Magellanic Cloud, we found a patchy struc- 
ture with the two major concentrations of late-type stars displaced 
a few (< 3) kiloparsecs further away. However, the vertical widths 
of red giant P-L relations seem to be significantly higher in the 
SMC, which can be attributed to the larger spatial depth range of 
the galaxy. Because of this, the calculated depth change is likely to 
be seriously affected and the interpretation of the results should be 
made with care. 
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